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Superconducting Fluctuation Effects on the Electron Spin Susceptibility in 

YBa 2 Cu 3 6 .95 
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The electronic spin susceptibility of YBa2Cu30e.95 has been measured with high precision up to 
24 Tesla with O nuclear magnetic resonance. Rs temperature dependence can be accounted for by 
superconducting fluctuations that result in a smooth crossover from the normal to the vortex liquid 
state. A magnetic field-temperature phase diagram for this crossover has been established having 
strong upward curvature. 



PACS numbers: 74.25.-q, 74.25.Nf, 74.72. Bk 

The upper critical field is large in high-T c materials be- 
cause of their small superconducting coherence lengths. 
However, a precise determination of H C 2(T) at high mag- 
netic fields is difficult because there is no well-defined 
signature of a phase transition. In magnetic fields sig- 
nificantly larger than the lower critical field (iJ c i(0) « 
100 G) the transition is broadened by the opening of a 
pseudogap in the electronic excitation spectrum. We per- 
formed a detailed experimental and theoretical study of 
the fluctuation effects on the spin susceptibility in op- 
timally doped YBCO and developed a quantitative un- 
derstanding for the onset of superconductivity solely in 
terms of a pairing pseudogap. This allows us to define a 
crossover field, iJ c2 (T), up to high magnetic fields. 

In high-T c materials the effects of a magnetic field on 
pairing fluctuations have beep discussed for the-diamag- 
netic response,D0 resistivitycl and heat capacity.!] In gen- 
eral, efforts to determine H C 2(T) have relied on ad hoc 
criteria that are not related to superconducting fluctua- 
tions. The specific heat has been analyzed by Roulin et 
al.a using several criteria to determine H C 2(T). SQUID- 
based magnetization measurements of Welp et aln show 
rounding in M(T) near the expected transition temper- 
ature as the field is increased. In this case the transi- 
tion temperature was determined by linear extrapolation 
of the temperature dependent diamagnetism in the su- 
perconducting state to intersect with the normal state 
magnetization. 

In the present work we determine H C 2(T) from mea- 
surements of the Pauli spin susceptibility, Xs- We mea- 
sure the 17 nuclear magnetic resonance (NMR) Knight 
shift and isolate the contribution of Xsi taking into ac- 
count orbital, diamagnetic, and vortex shifts. Because of 
the high precision of 17 NMR, we measure the temper- 
ature dependence of \s to better than 0.1% of the total 
normal state Xs value. Thus, we are able to make a quan- 
titative comparison of the data with the theory of pairing 
fluctuation corrections to Xs in the normal state. In the 
superconducting state there is a mean-field region where 



Xs is linear in temperature. A crossing-point analysis de- 
fines the crossover temperature, T C (H). We performed 
this analysis in magnetic fields up to 24 T. 

Our aligned powder sample is the same as that studied 
previouslyQ'cl having 30-40% 17 0-enriched YBa 2 Cu306.95 
prepared by solid state reaction. Low field magnetiza- 
tion data show a sharp T c at 92.5 K. The 17 NMR 
spectra were obtained from the fast Fourier transform 
of a Hahn echo sequence: 7r/2-r-7r-acquire and only the 
(1/2 «-> —1/2) transition of the 0(2,3) sites was stud- 
ied. High RF power allowed the use of short (rj 1.5 /is) 
7r/2 pulses (2.5 /xs at 2.1 T), giving a useful bandwidth 
> 100 kHz. The 17 0(2,3) (1/2 <-> -1/2) resonance has a 
low frequency tail owing to oxygen deficiency in a small 
portion of the sampled Its effect on our measurements 
was eliminated by performing a nonlinear least-squares 
fit in the frequency domain which isolates the dominant, 
narrow spectrum of optimally doped YBCO. Tempera- 
ture stability was ±0.1 K or better. For H a < 14.9 T 
we used superconducting magnets. The shifts measured 
at 8.4 T were the most precise. The high field measure- 
ments, 18.7 T to 24 T, were performed in a Bitter magnet 
adapted for modest homogeneity NMR at the National 
High Magnetic Field Laboratory in Tallahassee, Florida. 
High voltage arcing and temperature stability in the con- 
strained space of thq-Bitter magnet required that we de- 
sign a special probe.Q 

Excluding temperature-independent quadrupolar 
terms, the frequency of 17 NMR in YBCO can be 
written, 
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The gyromagnetic ratio of 17 is 17 7; H a is the static 
magnetic field applied parallel to the c-axis; AH V is a 
spatially varying field that results from pinned vortices, 
appearing below the melting or pinning temperaturefl 
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FIG. 1. The 17 0(2,3) spin shift at different magnetic 
fields for YBCO6.95. The spin shift is proportional to the 
Pauli spin susceptibility. Lines are guides to the eye. 

The orbital shift, 17 K or ^, is small and temperature inde- 
pendent. The diamagnetic contributions, AH^ & , whose 
origin is from surface supercurrents, are negligible at 
large fields B In the vortex liquid state the inhomoge- 
neous vortex fields, AH V , are motionally averaged. At 
low temperatures vortices are pinned, and the resonance 
for 17 nuclei in their vicinity is shifted to lower frequen- 
cies by 17 -fAH v compared to nuclei in the vortex liquid 
phase. This shift allows us to identify the 17 resonance 
associated with the vortex liquid B The temperature de- 
pendence of this part of the spectrum is given by 17 Kspini 
which is proportional to the electron spin susceptibility, 
Xs- The high sensitivity of 17 0(2,3) NMR to Xs, via the 
c-axis hyperfine coupling, is an advantage compared with 
that of 63 Cu(2) where the c-axis coupling is quite small 
and the copper resonance is broad. 

In Fig. [I] we show the temperature-dependent shifts 
for 17 0(2,3) measured over a magnetic field range of 2.1 
T to 24 T. The normal state values at 120 K are fixed to 
the value of 17 K S p- m in the normal state, 0.16% ± 0.01 %. 
17 Kg pm is the percentage spin shift relative to the Larmor 
frequency of the 17 nucleus and is attributable only to 
the electron spin susceptibility. The normal state value is 
determined independently from the extrapolation to high 
field of the difference 17 K(T = 100 K) - 17 K(T = 20 K). 
This expression extrapolates to the normal state 17 ^ S pin 
because other contributions to the total spectrum shift 
are either negligible at high fields (AH^ and AH V ), 
temperature independent ( 17 ^ rb)' or both (quadrupo- 
lar terms). The decrease of 17 K sp i n with decreasing 
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FIG. 2. 17 0(2,3) spin shift at 8.4 T for H\\c. The the- 
oretical calculation taking into account 2D pairing fluctu- 
ations is shown as curves, a) The inset shows a factor of 
7 expanded vertical scale demonstrating excellent agree- 
ment between theory and experiment. The crossed lines 
determine T c , as discussed in the text, b) The inverse of 
the derivative, (ci 17 A' S p m /dr)~ 1 , as discussed in the text. 
The dashed curves indicate calculations for T cm / = 72 K, 
75 K, 78 K, 80.9 K, 84 K, 87 K, 90 K, 92.5 K. 



temperature is smooth, showing no discontinuities in ci- 
ther magnitude or slope, suggestive of a crossover region. 
The crossover shifts to lower temperatures as the field is 
increased. 

Peak frequencies were determined by nonlinear least- 
squares fits to a gaussian around the peak region. In our 
stable superconducting magnet (8.4 T) this method al- 
lowed us to achieve precision in determination of the peak 
frequency of 1 part per million, corresponding to better 
than one part in a thousand precision relative to the to- 
tal normal state value of ifgpin- For data obtained with 
the Bitter magnet, corrections are required for variations 
of magnetic field with cooling water temperature, and 
thus the precision of the peak frequencies is only 5 parts 
per million. In Fig. Ha we plot 17 K spin at 8.4 T with an 



expanded version in the inset but on the same tempera- 
ture scale. The 17 Kgp m data have a monotonic decrease 
easily discernible below 110 K. 

The precision of our spin shift measurements at 8.4 T 
allows a quantitative comparison with the theory of su- 
perconducting fluctuations. Pairing fluctuation correc- 
tions to the Pauli spin susceptibility in the zero field 
limit tftve been considered theoretically by several au- 
thors En Gaussian fluctuation corrections diverge at the 
mean- field transition temperature T cm f. 
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FIG. 3. Feynman diagrams for the pairing fluctuation 
corrections to the Pauli spin susceptibility, to leading order 
in T cm f /Ef, as discussed in the text. 

Calculations for fluctuation contributions to the spin sus- 
ceptibility based on long wavelength static fluctuations at 
zero, field near T cm f predict (dSxs/dT)^ 1 oc T — T cm f in 
2Djj which cannot explain the curvature in our measurc- 
ments^hown in Fig. ||b, and (dSxs/dT)^ 1 oc */T — T cm j 
in 3D, El producing curvature opposite to that of our mea- 
surements. 

We calculate Xs in the weak-coupling limit for a quasi- 
2D, d-wave superconductor taking into account Landau 
quantization of the orbital motion of pairs by a mag- 
netic field (iJ||c). Our calculations include dynamical 
pairing fluctuations and we sum over all Landau levels 
in order to extend the range of validity of the theory to 
higher fields and temperatures. A detailed presentation 
of the theory is given in Eschrig et here we pro- 

vide a short summary. The pair fluctuation propagator 
for d-wave pairing with coupling constant g is given by 
HQ)' 1 = 9^-TJ2e n B 2 {tn,Q), where B 2 (e n ,Q) is an 
impurity renormalized two-particle susceptibility in the 
d-wave pairing chanpel as discussed in Mitrovic et aZ.L3 
and in Eschrig et a/Jll We use the notation Q = (u>i,q), 
which combines Matsubara energy, oji = 2nlT, and pair 
momentum, <f, of the fluctuation mode; the latter is quan- 
tized in a magnetic field. The results we obtain derive 
from the summation of all leading order pairing fluctu- 
ation corrections in T cm f / Ep, (Ep is the Fermi energy 
renormalized by exchange field corrections), given by the 
diagrams shown in Fig. [| (a-e), the Maki-Thompson (a), 
the density-of-states (b,c), and the Aslamazov-Larkin 
(d,e) contributions. K denotes the fluctuation propa- 
gator renormalized by external vertex corrections due to 
impurities. In contrast to the fluctuation corrections to 
the spin-lattice relaxation rate the contributions (d,e) in 
Fig. |3] have the same order in T cm f/Ep as (a-c). How- 
ever, they contain only one singlet pair fluctuation mode, 
K, the other mode in the particle-particle channel is a 
triplet impurity Cooperon-like mode, C. In the clean 
limit diagrams (d,e) vanish, in the dirty limit they give 
the main contribution. For intermediate impurity scat- 
tering all diagrams contribute significantly. It is possible 
to write the sum of all diagrams in Fig. in a compact 
way. Considering Bi(e n , Q) as a functional of the quasi- 
particle impurity self energy S(e„), the sum of all leading 



order pairing fluctuation corrections to Xs can be written 
as,rJ 

^CWE^^W). (1) 
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In Fig. H our calculation for 8.4 T is compared with 
experiment. The parameters extracted from the fit yield 
E F = 930 ± 30 meV, and T cm/ = 80.9 ± 0.3 K. We 
used the same scattering parameters as in our compari- 
son of the theory of dynamical pairing fluctuations with 
the field dependence of spin-lattice relaxation in Mitrovic 
et altB Dynamical fluctuations and orbital quantization 
effects produce the curvature shown in Fig. ^|b. The fit 
to our theory is performed in the region T > 90 K, and 
is shown by the heavy solid curve in Fig. |^a. Extension 
of the same fit to lower temperatures, as indicated by 
the thin solid curve, demonstrates that the theory fits 
the data well down to T = 85 K. Below this temperature 
critical fluctuations become significant. 

From the temperature dependence of the spin suscep- 
tibility in the mean field regime of the superconducting 
state it is possible to extrapolate linearly back to the sus- 
ceptibility of the normal state to determine a crossover 
temperature, T c . We determined the linear temperature 
dependence of the NMR spin shift in the superconduct- 
ing state from the maximum slope, d( K S p- m )/dT, and 
performed an extrapolation to the normal state shift of 
0.16% to find T c ; see Fig. |^a. The results are presented 
in Fig. |. In particular, at 8.4 T, we found T c = 84 K 
±0.5 K. The slopes we used are shown in the inset of 
Fig. | for various fields. Because the slope is expected to 
scale with A 2 , we can confirm the validity of our approach 
by comparing the slopes in Fig. ^ with direct measure- 
ments of the energ#-|gap, A. Tunneling and photoemis- 
sion measurements^ suggest a gap of 2A/fcs7c = 6 ± 1 
for YBCO. Using the mean field result for Xs without 
exchange corrections, but allowing for gap anisotropy, 
we estimate the scaled slope for d-wave pairing to be 
in the range 2.1 to 4.2, shown in the inset of Fig. [I| as 
a cross-hatched region. The anomalously large slope at 
2.1 T can be attributed to diamagnetic contributions to 
the frequency shift measurement. For all other fields the 
experimental slopes presented in the inset of Fig. || are 
consistent with our expectation for the slope of Xs- The 
crossover line, H C 2(T), in Fig. [| exhibits upward curva- 
ture. Qualitatively similar behavior, but at lower, fields, 
has been reported in specific heat experiments. □ With 
decreasing temperature there is a smooth crossover from 
the fluctuation regime to a vortex liquid phase. Devi- 
ations of our theory from experiment in Fig. |^ below 
85 K indicate onset of critical fluctuations. The smooth 
crossover suggests a relationship between critical fluctu- 
ations and the vortex liquid phase. 
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FIG. 4. H — T phase diagram from electron spin sus- 
ceptibility as determined by 17 NMR for YBC0 6 .9 5 . 
The dark squares represent a determination of H C 2(T) 
as discussed in the text. The curve is a guide to 
the eyei. Vortices become pinned in the cross-hatched 
region The open circle indicates onset of vortex 
pinning, an, extension of the earlier work of Bach- 
man et alu Dpen diamonds are from diamagnetism 
measurements £l Inset: Magnetic field dependence of 
(S = [d( 17 K spin )/dT] max x [T c / 17 ^ spin (120K)]). 

However, there is no adequate theory for the interplay 
between critical fluctuations and the fluctuating currents 
in the vortex liquid. 

The vortex liquid phase becomes well established at 
temperatures less than those given by H&iT). At lower 
temperatures vortices become pinned, as we have shown 
previouslyE Using two independent methods we deter- 
mined a region of the phase diagram in which pinned 
vortices are present, which we show as the hatched re- 
gion of Fig. l| We have extended that work to 27.3 T, 
shown as an open circle in Fig. |], where the onset of 
vortex pinning is observed by spin-spin relaxation. For 
low fields we point oui-that the pinning temperature and 
melting tempcratureou of untwinncd, single crystals co- 
incide. Our understanding of vortex pinning in combi- 
nation with our new measurements of H C 2(T) show that 
the region of liquid vortex matter is restricted to both 
high temperatures and low magnetic fields. 

In summary, we find that 17 NMR Knight shifts give 
a precise determination of the temperature dependence 
of the Pauli spin susceptibility, \s- We find significant 
rounding near T c indicating that superconducting fluctu- 
ations smear the transition. Consequently, the transition 
is best represented as a crossover from normal state be- 
havior to that of a vortex liquid. Our calculations for su- 
perconducting fluctuations, taking into account dynami- 
cal pairing fluctuations and the effects of orbital quanti- 
zation, are in excellent agreement with experiment in the 
temperature range down to T c . The decreasing suscepti- 
bility with decreasing temperature above T c can be fully 
accounted for by the opening of a pairing pseudogap. Fi- 



nally, we have established the H C 2(T) phase diagram for 
YBa 2 Cu 3 6 .95 up to 24 T. 
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